Effects of dilute Zn impurities on the uniform magnetic susceptibility of YBa 2 Cu 3 
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The effects of dilute Zn impurities on the uniform magnetic susceptibility are calculated in the 
normal metallic state for a model of the spin fluctuations of the layered cuprates. It is shown that 
scatterings from extended impurity potentials can lead to a coupling of the q ~ (n,n) and the 
q ~ components of the magnetic susceptibility x(q). Within the presence of antiferromagnetic 
correlations, this coupling can enhance the uniform susceptibility. The implications of this result 
for the experimental data on Zn substituted YBa2Cu307_i are discussed. 



PACS Numbers: 76.60.-k, 74.62. Dh, 75.40. Cx, 74.72. Bk 



Experiments on Zn substituted cuprates provide valu- 
able information on the magnetic properties of the host 
materials. Measurements have shown that Zn impurities 
cause an enhancement of the uniform magnetic suscepti- 
bility Qg] . From the nearly perfect Curie- Weiss temper- 
ature dependence of the enhancement, the magnitude of 
the effective moment, /i e ff, forming around the impurities 
has been extracted. In underdoped YBa2Cu306.66, Meff 
is close to jttBj while in YBa2Cu307 it is smaller by about 
a factor of 2.5. Neutron scattering experiments have also 
shown that Zn impurities modify the spectrum of the 
magnetic fluctuations in these systems ||||. Numerical 
calculations have emphasized the importance of the anti- 
ferromagnetic correlations of the host in determining the 
effects of Zn substitution H|. The effects of Zn impu- 
rities have been also considered in various spin-gapped 
and spin-gapless antiferromagnetic models H. Further- 
more, studies have been carried out for the spin-gapped 
underdoped phase of the cuprates fipf as well as for the 
d-w&ve superconducting phase ||,|lCj. Here results will 
be presented for the normal state of the layered cuprates 
where there are short-range antiferromagnetic correla- 
tions. Specifically, the effects of the impurity scatter- 
ing on the magnetic susceptibility will be calculated first 
for a single impurity and then it will be scaled to an 
impurity concentration of 0.5% in the dilute limit. It 
will be shown that scatterings from an extended impu- 
rity potential can lead to a coupling of the q ~ and the 
q ~ (tt, n) components of x(q)- Because of this coupling, 
the uniform susceptibility can get enhanced as antiferro- 
magnetic correlations develop in the system. In order to 
have a better understanding of this process, results will 
be presented on the related problem of how the q ~ 
and the q ~ (n, n) components of %(q) get coupled by a 
staggered charge-density-wave (CDW) field. 

The starting point is the two-dimensional Hubbard 
model with an additional term representing the inter- 
action of the electrons with a single impurity located at 
site r , 



H = -t £ (4 

(id) 



'^24<7 C iT + Yl ^ff(r ,ri)4cj ff . (1) 



Here Ci a (c ia ) annihilates (creates) an electron with 
spin a at site r^, t is the near-neighbor hopping ma- 
trix element, U is the onsite Coulomb repulsion, and 
fi is the chemical potential. The effective impurity- 
electron scattering potential is given by V c s (r , r) = 
Si/ ^ Sp„ H r > r o + Pv)> where p v denotes the sites at 
a distance v from the impurity site ro. Here, the effec- 
tive interaction is assumed to be static with a finite range 
extending a few lattice spacings away from the impurity. 
The importance of using an extended impurity potential 
has been previously noted |Il],ll2|]. For simplicity, in the 
following the hopping t and the lattice constant a will be 
set to 1. 

The single-particle Green's function is given by 

rP 
G(r i) T j ,iw n ) = - dre^ T (c ia (T)cU0)), (2) 



where u> n = (2n + 1)ttT. In the pure system with U = 0, 
one has in wavevector space Go(p,iuj n ) — (iui n — Ep)^ 1 
where e p = —2t(cosp x + cosp y ) — /z. If a single impurity 
is introduced at site ro, then 

G(r, r', iuj n ) = G (r, r', icu n ) 
+ Y, G (r,r",iuj n )T(r",r"' ) iu; n )G(r , ",r\iu) n ), (3) 

where the T-matrix for the impurity scattering is 

T(r,r',iuj n )=5(r,r')V e s{ro,r) 
+ J2 ^cff (r , r)G (r, r", w n )T(r", r', iu n ). (4) 

r" 

This is illustrated diagrammatically in Fig. p](a). The 
magnetic susceptibility is defined as 

X (r,r',mv)=/ dr e ^ T (m-(r,r)m+(r',0)), (5) 
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+ cl J(J c la ) + u}_^ c,; T c iT cj iCli 



where m + (r) = cLc r |, m (r) = cl^c T ^, and Lu m = 
2rmrT. First, the effects of a single impurity will be 



calculated for U — 0, giving the irreducible suscepti- 
bility xo, and then the effects of the Coulomb correla- 
tions will be included. The diagrams representing the 
effects of a single impurity are shown in Figs. 0(b) and 
(c). Both the self-energy and the vertex corrections 
need to be included [J13), and the resulting expression 
for xo(r,r') = Xo(r,r',«w m = 0) is 



(r,r') = -T^2 \ [G(r,r / ,tw n )Go(r / ,r,tw n ) 



+G (r, r',iu) n )G(r',T,iu> n ) - (G (r,r',iw„)) 2 ] 
- ^2 G (r,ri,iw„)Go(r2,r',«w n )T(ri,r 2 ,w n ) 

ri,r 2 ,r 3 ,r 4 

xr(r3,r4,iw„)Go(r3,r,iw„)Go(r',r4,iw n ) }. (6) 



Here the self-energy corrections are included by summing 
over the terms in the square brackets rather than simply 
summing over GG. In addition, the irreducible impurity- 
scattering vertex has been used in calculating the vertex 
corrections to xo instead of the reducible one. These 
are necessary in order to prevent double counting when 
calculating the effects of just one impurity on xo- 

The Coulomb correlations are included by using the 
random-phase approximation, 

X (r, r') = X o(r, r') + UV] X o(r, r")x(r", r'). (7) 



Upon solving for x(r, r')i the Fourier transform is taken, 
leading to xlQil')- At this point the impurity averag- 
ing can be done. This restores the translational invari- 
ance, and the diagonal susceptibility x(q) = ^qq'X^q') 
is obtained. In order to minimize the finite size effects 
on x(q)j t ne single-particle Green's functions used in 
Eqs. (||) and (||) are evaluated on large lattices basically 
without any finite size effects. Equation (Q), on the other 
hand, is solved on smaller L x L lattices for a single im- 
purity located at the center, but the resulting x(q) i s 



(a) 



(b) 
I. 



(c) 



* 



i i \ 



■ + c ' c + c ' < » ( + c '< ' <:' < + ■ 



ooo<x> 



<s> - Ci>a><2>- 

FIG. 1. Feynman diagrams for (a) the dressed sin- 

gle-particle Green's function G, and (b)-(c) the irreducible 
susceptibility xo within the presence of a single impurity. 



scaled to the L — > oo limit while keeping the impurity 
concentration, ni mp , fixed. This is an efficient way of con- 
trolling the finite size effects when working in the dilute 
limit. The following results on x(q) are for n; mp = 0.005, 
U = 2t, electron filling (n) = 0.86, and T = 0.1. For 
these values of U and (n), the pure system has short- 
range antiferromagnetic fluctuations. In addition, the 
impurity potential was assumed to have a range of 2 lat- 
tice spacings with the following parameters: Vo = —20, 
Vi = 0.5, V^ = -0.5, and V 2 = -0.25 @. 

Figure 0(a) compares x(l) versus q for the pure and 
the 0.5% Zn substituted cases obtained by using a 14 x 14 
lattice in Eq. (0). Figure [|(b) shows results on x(q) ver- 
sus q x along q x = q y . Here the solid line represents x(q) 
for an infinite pure lattice, and the open circles are for 
the 14 x 14 pure lattice. One sees that the finite size ef- 
fects are small. The filled circles are the results obtained 



on the 14 x 14 lattice for 



nmp 



= 0.005. Comparing with 



the open circles, one observes that x(q —> 0) is enhanced 
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FIG. 2. (a) x(q) versus q along the path shown in the in- 
set for the pure and the 0.5% Zn substituted systems on the 
14 x 14 lattice, (b) x(c[) versus q x along q x = q y . The solid 
line is for the pure case on an infinite lattice. The open and 
the filled circles are for the pure and the impure cases, respec- 
tively, on the 14 x 14 lattice. One observes that x(q ~^ 0) is 
enhanced by the impurity scattering. The triangles are also 
for the impure case on the 14 x 14 lattice, but now the im- 
purity averaging has been done before including the Coulomb 
correlations. This neglects the mixing of the q ~ and the 
q ~ (n, n) components of x(q)> an d subsequently x(q) is sup- 
pressed. 



by the impurity scattering. This can be understood bet- 
ter by considering Eq. (ffl) in wavevector space, 

x(q, q') = xo(q, q') + tf £ xo(q, q")x(q", q')- (8) 

q" 

Here, one sees that once the translational invariancc 
is lost, the Coulomb interaction mixes the different 
wavevector components of the magnetic susceptibility. 
This allows for the q ~ component to be influenced by 
the antiferromagnetic correlations. Through this mixing, 
x(Q = (7T,7r)) also gets enhanced as seen in Fig. ||(a). 
The open triangles in Fig. ||(b) are for the impure case, 
too, but here the impurity averaging has been done be- 
fore including the Coulomb correlations. This does not 
allow for the mixing of the q ~ and the q ~ (x, tt) 
components, and consequently x(q) is suppressed. 

Taking the difference of the filled and the empty cir- 
cles shown in Fig. 0(b), one obtains the enhancement 
of the magnetic susceptibility by the impurity scatter- 
ing, Ax(q). In Fig. ||(a), A%(q) versus q x is shown as 
the system size is scaled from 14 x 14 to 28 x 28 while 
keeping roi mp fixed at 0.005. One observes that the finite 
size effects on Ax(q) are negligible. The inverse of the 
peak width of A%(q) gives the size of the region around 
the impurity which contributes to the enhancement of 
the uniform susceptibility. For comparison, A%(q) ob- 
tained using an onsite impurity potential of zero range 
with Vo = — 20 is shown by crosses in Fig. @(a). In this 
case, A%(q) is significantly smaller and it is featureless 
in q. Figure 0(b) shows the temperature evolution of 
Ax(q) for the extended potential. Here the growth of 
^x(q — > 0) with decreasing T is seen. On the other 
hand, for the onsite potential Ax(q) has a weak T de- 
pendence (not shown here). 

In order to display the coupling of the q ~ and 
the q ~ (x, n) components of x(q) by impurity scat- 
tering, the Coulomb repulsion U entering Eq. (|7j) has 
been varied by small amounts while keeping the rest of 
the parameters fixed. The resulting Ax(q — > 0) ver- 
sus x(Q = (^i 7 !")) is plotted in Fig. [|(c), where a lin- 
ear dependence between these two quantities is seen. 
For instance, by increasing U by 4% from 2.0 to 2.08, 
x(Q = (tv, tt)) increases from 3.9 to 5.7, while Ax(q — ► 0) 
doubles. These calculations were repeated using various 
other sets of V v 's and similar results were obtained for 
the coupling between the q ~ and the q ~ (tt, tt) com- 
ponents of x(q)- 

The enhancement of x(q ~* 0) seen m Fig- PK a ) is 
about 2%. However, a direct comparison with the exper- 
imental data was not carried out because of the simplicity 
of the model. For instance, the exact values of V„'s are 
not known, and the effects of the Coulomb correlations 
on the single-particle Green's functions and the T de- 
pendence of V e s are not taken into account. For these 
reasons, it is not possible to make a direct comparison 
with the Curie- Weiss behavior, either. 
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FIG. 3. (a) Finite size effects on Ax(q) versus q x along 
q x — q y for the extended impurity potential. For comparison, 
Ax(q) obtained using an onsite impurity potential is shown 
by crosses, in which case Ax(q — » 0) is smaller, (b) Tem- 
perature evolution of Ax(q) versus q x . (c) Ax(q — > 0) versus 
x(Q = (7r, Ti")) obtained by varying U by small amounts, which 
shows that these two quantities are correlated. 

These results on the effects of scattering from ex- 
tended impurity potentials can be understood better if 
one considers how a staggered CDW field affects x(q) 
in a pure system. This is a simple example which in- 
volves the mixing of the q and the q + (tt, tt) compo- 
nents of x(q)- I n fact, the calculation of x(q) for a single 
impurity reduces to this problem, if one only keeps the 
component of the effective impurity-electron interaction 
which transfers Q = (tt, tt) to the quasiparticles. Here, 
one begins by introducing a staggered field A coupling 
to the site occupation number, ^2 icr Ac^c^ exp (iQ • r,) 
with Q = (tt, tt). The resulting irreducible susceptibil- 
ity has both diagonal and off-diagonal terms, Xo(q, q) 
and Xo(q + Qjq)> which are shown diagrammatically 
in Fig. 0(a). Within the random-phase approximation, 
Xo(q, q) and xo(q + Q, q + Q) are coupled through the 
off-diagonal term xo(q + Q, q) leading for x(q, q) to 



X(q,q) = J D- 1 {xo(q,q)(l 



C/Xo(q + Q,q + Q)) 

+^xo(q + Q,q)}, 



(9) 



where D = (1 - C/ X o(q, q))(l - £/Xo(q + Q, q + Q)) - 
t^ 2 Xo(q + Q>q)- This calculation of x within a CDW 
field is similar to that within a spin-density-wave field 
Jl5[ . The irreducible susceptibilities seen in Fig. |j(a) 
can be evaluated and the resulting x(q>q) is displayed 
in Fig. 0(b), where the enhancement of x(q, q) by A is 
seen. It can be shown that, for small A, the off-diagonal 
term Xo(q + Q: q) remains finite in the limit q — > [Gil, 
and it is approximately equal to (A//i)xo(q, q), while 
Xo(q, q) is changed little by A. Hence, it is expected from 
Eq. (pi) that the uniform susceptibility will get enhanced 
by the antiferromagnetic correlations, as A is turned on. 
In addition, it can be shown that Xo(q+Qi q) is a smooth 
function of q for q ~ 0. Hence, from Eq. (||) one observes 
that the q structure of the enhancement in x(q, q) for 
q ~ reflects the structure of x(q^q) f° r q ~ (^j 71 ")) 
which is slightly incommensurate at this temperature. 

It is important to note that here the enhancement of 
the uniform susceptibility is not due to an enhancement 
of Xo(q> q) by A, since these small values of A have little 
effect on Xo(q, q)- Rather, it is due to the nonvanishing 
of the anomalous susceptibility Xo(q+Q> q), which allows 
for a coupling to the antiferromagnetic correlations. In 
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the impurity problem, the enhancement of x(q ~ 0) 
occurs similarly. In Eq. (||), the off-diagonal terms 
Xo(q, q" 7^ q), which are nonvanishing since the impu- 
rity scattering does not conserve momentum, couple the 
q ~ and the q" ~ (71", 7r) components of x, enhanc- 
ing the uniform susceptibility. This coupling is strong 
especially when an extended impurity potential is used. 
However, if the impurity averaging is done before includ- 
ing the Coulomb correlations, this coupling is neglected 
and x(q) is suppressed with respect to the pure case, as 
it was seen in Fig. 0(b). 

In summary, the effects of dilute Zn impurities on the 
uniform magnetic susceptibility have been calculated in a 
metallic model which has short-range antiferromagnetic 
correlations, but does not necessarily have a spin gap. 
It has been found that impurity scattering through an 
extended potential leads to the mixing of the q ~ (tt, it) 
and the q ~ components of x(q). Because of this cou- 
pling, the antiferromagnetic correlations can enhance the 
uniform susceptibility. The microscopic model presented 
here predicts a strong dependence of Ax(q — * 0) on 
the strength of the antiferromagnetic correlations. This 
could play a role in determining the hole-doping and the 
temperature dependence of Ax(q — ► 0) in Zn substituted 
YBa2Cu307_,5. However, it must be kept in mind that 
these results depend on the nature of the effective impu- 
rity interaction and on the validity of the approach used 
to treat the Coulomb correlations. 

The author thanks H.F. Fong for helpful discussions. 
The numerical computations reported in this paper were 
performed at the Center for Information Technology at 
Kog University. 
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FIG. 4. (a) Feynman diagrams representing the diagonal 
and the off-diagonal irreducible susceptibilities xo(q, q) an d 
Xo(q+ Q, q), respectively, within the presence of a staggered 
CDW field, (b) x((\, q) versus q x along q x — q y for various 
values of the CDW amplitude A, where the enhancement of 
X by A is seen. 
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